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1.  INTRODUCTION 

This  report  describes  a  software  package  for  the  generation  of  range  height-angle  charts, 
which  are  commonly  used  to  provide  qualitative  guidance  on  the  propagation  characteristics 
of  electromagnetic  waves  through  the  atmosphere  in  the  solution  of  radar  and  communication 
problems.  Such  charts  show  the  interrelationship  between  range  to  a  target,  the  height  of 
the  target,  and  the  angle  that  a  ray  has  to  be  launched  from  the  source  to  reach  the  target 

position.  They  may  appear  in  many  forms,  but  most  often  have,  and  are  most  useful  when, 

the  range  and  height  axes  are  transformed  sc  that  the  propagation  paths  are  straight  lines, 

Tlie  need  for  such  charts  arises  because  within  the  atmosphere  the  paths  are  governed  by  the 

spatial  variation  of  the  atmospheric  refractive  index,  and  are  generally  curved  downwards 
due  to  the  normal  reduction  in  refractive  index  with  height.  Under  more  extreme 
circumstances  this  curvature  can  match  the  curvature  of  the  Earth,  and  consequently  the 
radar  waves  can  propagate  considerably  beyond  the  geometric  horizon,  a  condition  known  as 
ducting.  Thus,  most  often,  given  a  specified  target  position  and  a  radar  position,  it  is 

not  immediately  apparent  if  a  propagation  path  exists  between  them  for  a  particular 
atmospheric  refractive  index  profile.  Nor  is  it  simple  to  deduce  the  launch  angle  of  the 
ray  from  the  radar  that  passes  through  the  t<-  3d  position,  and  it  is  in  answer  to  such 

questions  that  range-height-angle  charts  are  commonly  constructed. 

Most  charts  that  have  been  published  to  datefref.  1 ,2)  usually  assume  a  standard  atmospheric 
refractive  index  profile(ref.5),  and  place  the  source  at  ground  level.  Although  it  is 
possible  to  deduce  if  a  ray  path  exists  between  two  elevated  points  by  determining 

separately  the  horizon  distance  of  each  of  the  points  and  comparing  the  sum  of  these  with 
the  physical  separation  of  the  points,  the  launch  and  arrival  angles  of  the  ray  between  the 
points  cannot  be  determined  in  this  process.  The  alternative  of  providing  an  extensive  set 
of  charts  covering  the  range  of  source  heights  of  interest,  typically  from  0  to  100, (XX)  ft 
for  most  aircraft  applications,  and  for  various  refractive  index  profiles  is  not 

particularly  attractive,  but  with  the  advent  of  high  speed  digital  computers  with  good 
graphics  capabilities,  it  is  possible  to  generate  charts  simply  for  each  specific  situation 
of  interest.  This  report  dv.ails  a  software  package  that  has  been  developed  for  this 
purpose  by  adapting  a  general  purpose  program  for  tracing  rays  through  the  atmoxphere(ref.4) 

that  was  written  for  use  with  a  radar  system  modelling  package  under  development  in 
Microwave  Radar  Division.  Reference  3  details  in  alternative  method  for  the  generation  of 
charts  for  elevated  sources,  with  a  standard  exponential  reference  atmosphere.  The  charts 
presented  therein  have  the  angular  ordinate  scaled  by  a  factor  of  10,  giving  them  an  unusual 
lorm  and  restricting  the  coverage  of  high  and  low  aspect  angles.  No  such  limitations  exist 

with  our  product. 
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The  software  program  that  has  been  developed  is  completely  general,  in  that  the  user 
specifics  the  source  height,  the  limits  in  range  and  height  to  which  the  rays  are  to  be 
traced,  and  the  functional  form  of  the  atmospheric  refractive  index  profile.  The  output  of 
the  program  is  a  range-height -angle  chart,  in  which,  by  appropriate  transformations  of  the 
range  and  height  axes,  the  ray  path  from  the  radar  is  depicted  as  a  straight  line. 
Previously  published  charts  have  commonly  depicted  ray  path  length  as  the  range  parameter, 
resulting  in  portions  of  ellipses  for  the  constant  iange  curves  at  short  range.  We  have 

however  chosen  to  use  ground  range  as  the  range  parameter,  resulting  in  a  slightly  different 
form  for  the  constant  range  curves  as  we  have  found  this  to  be  more  useful  for  our  purposes. 
At  moderate  range  the  differences  are  only  slight,  but  for  specific  requirements  where  the 
differences  are  significant,  it  would  be  straightforward  to  adapt  the  software  to  use  the 
ray  path  length  as  the  range  parameter. 


2.  THEORY 


2.t  Ray  tracing  through  the  atmosphere 

The  fundamental  requirement  in  the  generation  of  range  height  angle  charts  is  to  be  able  to 
trace  the  path  of  rays  as  they  propagate  through  the  atmosphere.  In  this  work  we  will  rely 
upon  the  techniques  of  geometrical  optics  for  this  task,  noting  as  an  aside  that  such 
techniques  will  be  inadequate  at  and  near  caustics  of  the  rays,  but  also  observing  that  this 
package  is  not  intended  to  be  used  for  refractive  index  profiles  which  give  rise  to  such 
phenomena.  The  theoretical  basis  for  ray  tracing  is  treated  more  fully  elscwhere(rcf,4), 
but  may  be  summarized  for  our  purposes  as  follows. 

(1)  At  any  point  in  the  (horizontally  stratified)  atmosphere,  the  local  curvature,  K,  of 

a  ray  will  be 


K 


1 

n 


cos  a 


dn 

dz 


where  a  is  the  angle  of  the  ray  to  the  horizontal,  and  n  is  the  atmospheric  refractive 
index. 


(2)  Assuming  that  the  curvature  of  the  ray  is  constant  along  a  sufficiently  small 
segment  of  the  ray  path  (say  100  m  for  example),  then  the  coordinates  of  the  end  point  of 
the  segment  may  be  simply  determined  from  the  coordinates  of  the  start  of  the  segment 
(see  figure  1):  thus  compute 
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fa)  Jr  =  A.JS 

(b) 

fc) 

(ci) 

(e) 

(0 

where  (h  ,  8  )  defines  the  height  above  the  surface  of  die  Earth  and  the  angular  position 

th 

based  on  the  centre  of  the  Earth  of  the  start  of  the  i  segment,  and  (h  .,8  )  the 

t  1  + 1  i  + 1 

same  parameters  at  the  end  of  the  i  segment,  ie  the  start  of  the  i  +  1  segment.  The 

angle  a  in  these  expressions  is  the  angle  of  the  ray  to  the  local  horizontal  to  the 

Earth’s  surface,  and  R  is  the  radius  of  the  Earth.  Figure  !  illustrates  the  relationship 

between  these  parameters. 

(3)  Using  the  above  procedure  iteratively,  rays  can  be  launched  from  the  source  and 
traced  until  either  the  range  or  height  exceeds  the  boundaries  specified  for  the  chart, 

or  the  ray  intersects  the  surface  of  the  Earth. 

2.2  Intersection  of  a  ray  with  a  constant  height  curve 

To  construct  the  constant  height  curves  of  the  chart,  we  need  to  determine  the  points  on 

each  of  the  rays  at  which  the  heights  of  the  curves  are  crossed.  The  problem  is  very 

similar  to  that  of  tracing  the  rays  outlined  in  the  previous  section,  except  that  in  this 

case,  we  know  the  final  height  of  the  end  point  of  the  segment  but  not  the  segment  length. 

Assuming  that  the  i'h  ray  path  segment  straddles  the  height,  H  for  which  a  curve  is 

required,  then  on  referring  to  figure  1,  in  which  h  is  replaced  by  H.  and  putting 

Z  =  h  +  R 

i  i 

Z  =  H  +  R 


.11,  = 


a  =  y  (rr  -  Jr) 


h  +  R  =  f(h  +  R)2  +  JL2  -  2(h  +  R)  JL  cos  (a  +  a  ) 
i+1  I  i  i  il 


AO  arcsin  i sin  («  +  <7  )  ,  . > 

i  ^  i  i  h  +  aj 


a  =  a  +  A 8  -  At 

i  +  $  i  i 


it  can  be  shown  that 
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7~  -  Z‘  +  AL“  +  2ZJ L  sin|«  -  — ^j- 

On  expanding  the  trigonometric  function,  and  noting  that 

2  sin  (Jr/2) 

K 

where  A'  is  the  curvature  of  the  ray,  a  quadratic  in  A L~  wliich  can  easily  be  solved,  is 

obtained.  The  form  is 

A  JL4  +  B  AL2  -  C  =  0 

in  which 

A  -  (A'Z)2  -  2KZ  +  1 
1  1 

B  =  2  ((Z2  -  Z2)(A'Z  cos  <v  -  1)  -  2Z2  sin2u  1 

l  .  i  .  ‘  ') 

and 

c  =  [z2  -  z2}2 

There  are  two  positive  real  solutions  to  this  equation,  of  which  the  one  we  seek  can  be 
shown  by  geometrical  construction  to  be  the  smaller.  The  pioblcm  of  whether  a  real  solution 

exists  does  not  arise  since  the  procedure  is  invoked  only  for  segments  which  straddle  the 
specified  height  curves.  After  solving  for  AL,  the  arc  iength  JS  can  be  determined,  and  from 
this,  the  ray  path  length  to,  end  the  coordinates  of,  the  intersection  point  on  the  chart 

may  be  calculated. 

2.3  Intersection  of  a  ray  with  a  constant  range  cun 

The  intersection  point  of  a  ray  with  a  particular  range  curve  follows  in  a  similar  fashion. 
Referring  again  to  figure  1,  and  assuming  that  the  i'h  segment  of  the  ray  path  crosses  the 
range  curve,  then  since  we  know  A 6.  from  the  range  and  since 

«  i  f  dL 

,1r  1  2  arc  sin  i  A 


SR  L -0027  -TM 


on  applying  the  sine  law  for  the  ratio  of  sides  of  a  triangle,  wc  obtain  the  result 

J  L  7  i__ 

sin  (  A  8  )  sin  y 

Z 

i 

cos  (AS  +  a  ■  i7  2) 

i  l 


i 

a  quadratic  for  JL~  similar  to  the  previous  case  can  then  be  deduced  in  which  the  parameters 
A,  B,  and  C  assume  the  values 


A 


1 

4 


B  =  -  Z  sin  8  sin(S  +  nt.)  K  -  cos' (t >  +  or.) 

i  i  i  1  ll 

C  =  Z2  sinV 

i  l 

Once  AL  is  determined,  At  follows,  and  from  it,  the  path  length  along  the  arc  to  the 

intersection  point  of  the  ray  path  with  the  specific  range  curve,  and  the  height  of  the 

intersection  point  may  he  calculated. 

2.4  Construction  of  range-height  angle  charts 

A  range  height  chart  may  be  constructed  by  tracing  the  paths  of  a  set  of  rays  launched  at 

various  angles  from  the  source,  and  generating  tables  of  the  path  lengths  of  the  rays  to  the 

intersection  points  of  the  rays  with  the  heights  and  ranges  for  which  curves  are  required. 

The  loci  of  the  points  of  constant  height  and  constant  generated  by  the  set  of  rays  provide 

a  variable  scale  for  the  axes,  whereby  the  ray  paths  may  be  drawn  as  straight  lines.  The 

procedure  foi  the  generation  of  the  range  and  height  contour  intersection  points  for  each 

ray  is  terminated  when  the  end  point  of  a  ray  segment  exceeds  either  the  maximum  height  or 

the  maximum  range  specification  or  intersects  the  surface  of  the  Earth.  To  provide  an 

accurate  boundary  for  ihe  chart,  the  end  point  of  the  last  segment  is  adjusted  to  terminate 

precisely  on  the  height  oi  range  boundary,  using  a  variation  of  the  procedures  described  in 

Sections  2.3  and  2.4  in  which  all  the  parameters  of  the  end  point  of  the  increment  are 

revised  once  the  step  increment  Ji,  is  known.  For  each  intersection  point,  the  mapping  used 

to  determine  the  location  on  the  chart  (X  ,  Y  1  with  respect  to  the  chart  origin  is; 

P  P 

X  -  R  cos  n  /  X 

1  p  -S 

Y  (R  Mil  (t  4  H  )  /  Y 

p  p  s  S 

where 

R  is 

p 


the  toial  ray  path  length  to  ihe  intersection  point 
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rt  is  the  launch  angle  of  the  ray  from  the  source 
H  is  the  height  of  the  source  above  the  surface  of  the  Earth 

X  ,  Y  are  the  x  and  y  axis  scale  parameters 

Determination  of  the  scale  factors  will  be  discussed  later. 

Given  then,  the  points  on  the  chart  at  which  the  rays  intersec*  each  of  the  specified 

heights  and  ranges,  the  curves  for  constant  height  and  range  can  be  constructed  by 

connecting  the  respective  points.  Since  the  rays  are  traced  at  various  starting  angles  from 

the  source,  beginning  at  +90°  and  extending  through  O'  to  -90°,  it  does  not  follow  that  the 
points  of  the  height  and  range  curves  are  generated  in  sequential  order:  in  particulai  for 

an  elevated  source  height  it  is  possible  for  rays  launched  below  ()'  to  cross  particular 

height  boundaries  twice  provided  that  they  do  noi  first  intersect  the  surface  of  the  Earth 
Connection  of  the  points  of  each  curve  in  the  order  of  their  generation  would  not  provide  a 
monotonic  result.  This  problem  has  been  resolved  in  the  software  by  sorting  the  points  foi 
each  height  curve  and  for  each  range  curve,  into  increasing  range  and  height  order 

respectively  before  plotting. 


The  upper  boundary 

of  the  chart  is  formed 

by  the  rays  which 

lermi"afe 

,i( 

the  height  1 

\y.ir. 

the  right  hand  side 

by  the  rays  that  terminate  at  the  larige 

limit,  and 

the 

bottom  be 

the 

rays  which  terminate 

on  the  surface  of  the 

Earth  and  by  the 

extension 

of 

the  tangent 

ray 

Certain  pathological 

conditions,  for  example 

only  having  one 

ray  tern. mating  on  the 

right 

hand  boundary,  are 

taken  care  of  by  making 

assumptions  -oout 

the  shape 

he  ;vr  v 

greater  detail  is  rot  available.  In  the  quoted  example,  the  right  hand  boundary  is  assumed 
to  be  vertical.  Where  there  are  no  rays  terminating  on  the  maximum  range  boundary,  the 

tangent  ray  and  the  maximum  height  boundary  intersect  each  other  Another  problem  dial 

arises  is  that  most  often  me  {Hunts  that  define  the  intersection  of  each  ot  the  range 
curves  with  the  maximum  height  boundary  (and  the  corr  pkmetitary  situation)  wdl  not  tv 
generated  in  the  ray  tracing  this  problem  has  been  rest  Ived  by  linearly  interpolating  the 
respective  curve  da‘a  to  provide  the  intersection  point  with  the  boundary  Provided  that 
the  successive  rays  are  closely  spaced  (in  chart  coordinates),  tins  gives  an  acceptably 
accurate  result 

A  major  difficulty  in  developing  the  software  was  experienced  in  getting  the  rays  suitably 

spaced  to  provide  an  adequate  number  of  data  (tonus  foi  each  range  and  height  curve  for  the 

cx|>ectc(i  practical  values  of  the  range  and  height  limits  'Hus  problem  has  been  resolved  by 
determining  approximately  the  launch  angles  Horn  the  source,  of  the  ray  dial  pusses  through 
tlic  range  height  boundary  intersection,  ac.J  the  ra,  tangent  to  the  Earth'  ,  suifacc.  using  a 
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4/3 rJ  Earth  radius  inode!.  From  the  geometry  of  figure  2,  it  follows  that  the  launch  angle  p 
;o  the  range-height  boundary  intersection  is: 


P 


tan 


+  H  )  cos  S  -  H 
mix  s 

( R  +  H  )  s  i  n  d 

c  f  f  in  t  x 


where 

R 

c  inn  x 


R  =  maximum  range  of  ray  trac  ig 

max 

R  =  effective  Earth’s  radius  —  4/3  times  actual  Earth’s  radius 
eff 

H  =  source  height 

U  ~  maximum  height  of  rav  tracing 

max 


Similariy  the  launch  angle  //  for  the  ray  tangent  to  the  Earth’s  surface  is  (from  figure  2) 


/>' 


c  f  t 
+  R 


ami  the  angle  y  to  the  bottom  corner  of  the  chart  is 

H  4-  (  1  -  cos  i5 /  R 

s  eft 

R  sin  S 

<•  if 

The  software  package  then  traces  ravs  n  various  intervals  t S ' ' .  10°  or  30°)  from  >  v  to  1()°. 

then  m  0  V  intervals  to  />.  0  1  intervals  to  O',  0  03'’  intervals  to  //,  i)  l'  n  rivals  to 
}\  1  intervals  to  10’,  and  various  intervals  to  In  fact  this  software  ,  slightly 

more  complicated  than  suggested  above  based  on  the  relative  values  of  the  break  points  for 
different  angular  steps,  and  the  values  of  //,  d  and  y  (sec  Listing  2).  However  the  general 
principles  arc  adequately  illustrated,  with  the  finest  division  of  launch  angle  being  used 

between  0  and  the  ray  tangent  to  the  surface.  In  general  tin  scheme  has  worked  well,  but 
for  jHissibly  extreme  uses  of  this  sob  ware  m  which  the  specified  maximum  height  is  very 

small  with  respect  to  the  maximum  range,  the  resulting  charts  may  not  have  enough  (Hunts  to 
give  an  adequately  .’"Moth  set  of  range  and  height  curves.  In  such  circumstances  it  will  lv 

necessary  to  modify  the  angular  imrcments  and  break  points  to  suit  the  requiiement  at  hand. 
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The  selection  of  the  X  and  Y  axis  scale  parameters  and  the  position  of  the  origin  so  that 

the  chart  occupies  the  entire  output  page  is  not  straightforward.  For  charts  with  extended 

ranges  (see  for  example  figure  6)  a  considerable  portion  of  the  chart  lies  below  the  origin, 

and  consequently  the  origin  has  to  be  raised  (relative  to  the  bottom  of  the  chart)  and  the  Y 

scale  compressed  to  fit  the  chart  in.  In  this  package  these  parameters  are  determined  using 

a  4/3 r<1  Earth  approximation  to  produce  a  chart  8”  x  6'  in  dimensions,  and  then  applied  to 

the  ray  tracing  output.  Two  situations  are  of  significance  in  these  calculations,  whether 

the  tangent  ray  to  the  surface  of  the  Earth  touches  the  surface  before  or  after  the  maximum 
range. 

If  the  maximum  range  is  less  than  the  tangent  range,  then,  referring  to  figure  2(a).  we 

compute. 


<t>  = 


from  which  the  scale  factors  X  .  Y  follow  immediately,  thus 

s  s 


X  R  sin  <i>  /  X 
s  rtf  I. 


Y  -In  +  R  (1 

s  I  man  rff 


COS  <£)>  /  Y 


where  X ^  and  Y  are  the  X  and  Y  dimensions  of  the  chart. 


If  the  maximum  range  is  greater  than  the  tangent  range,  then  determining  the  scale  factors 
is  a  little  more  complicated  Refemn;,  to  figure  Shi.  we  compute 


* 


■  i  I 


K  R 

mas  i 


I  he  height  of  the  tangent  rav  above  the  surface  of  the  Earth  at  range  K  is 

llllk 


H 

i 


R 

Ml 


I'll  S 

i  n  s 


i!  this  height  is  kircatci  than  H  then  the  an^le  ^  is  revised  so  that  the  undent  i.n, 

max 

k  riuiihili  s  ;u  »n.»x iinuiii  height,  thus 
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V  ~  cos 


R  +  H 

c f  f  in  ■ x 


The  scale  factors  for  the  axes  thus  become: 


■  R  -f  tanu/)  cos  <£  X 

s  rff  r 


and 


^  >  H  H  +  R  (tan  <t  f  t;i 


irmx  s  df 


:tn  y ;  sin  >  Y 


fhe  length  of  the  vertical  axis  then  follows: 


L  -  H  Y 

tn«x  i 


am!  from  this  .he  ros.non  of  ,he  origin  may  he  deterrmnod  .so  .hat  fhe  eh, an  f„s  neatly  on 
the  output  page.  In  our  .system.  <*e  use 


X  I  S 

o 


v  i. 


sshieli  roughly  centres  the  ,  hart  on  an  A  4  page 


3.  INVOKING  ANI)  TSINli  THK  SOFTWARE 


'  I  t\er  interface 


n,C  l,,'V  ""s  K'K,n  ^  »•»«'  "<  •<  '<rger  package  which  eau  carry  out  , 

':,rK,V  ,:,V  . . .  iHM  i?"  '>'-»«  n,e  . . ,w;1rc  lor  the  package 

contained  a,  a  ...nets  of  ,1  ,ta  sets  hcg.nn.ng  with  the  descriptors  Jl  W  RAYTRAOM  ’.  (he 

ll'C  '‘‘V  'nit  mg  ami  chart  geneiation  procedures  is  in 
M  V,  KAN  IKAl-l-  i'l  It  M  AIN  A  and  other  memlvrsol  JI  W  RAY  IKAti:  i'l  V  are  called  hv  MAIM  |1u. 

im.rlate  is  via  S  P  l-  panels,  invoked  initially  hv  the  I'l. 1ST.  Jl  W  II  IS  ft  R  A  YTR  AIT  )' 
H"'  ^•IS,‘  ,he  SPY  menus  ami  panels  Much  the  ,.,v  „,g 

-htsare.  ami  then  displays  the  l„s,  panel  to  select  the  r.,v  ttaeng  task  te, pored  II, e 
•  ’ptnui.s  available  ,n  this  package  at  preset  arc 


I  1  tr.ii  c  rax  \  in 


m  an  am 


•nna  f. •?  a 'i  mhMrarx  rcl/acuvc  index  pi,,filc 


"''c  ",e  ru  l‘cl'WC11  l-in.s  1  smt.it'le  only  tor  ,c|,.„mc  ,,,dcx 

ptolile  will,  h  docs  ill,,  ci'r  use  (.,  dll,  tihg' 
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( 3 ;  plot  range-height-anglc  charts 

(4)  plot  propagation  loss  between  two  points  as  a  function  of  height  of  the  receiver 

(5)  generate  resolution  cells  for  time  difference  height  finding  for  a  4/3  or 

other  effective  radius  atmospheric  profile 

In  the  case  of  the  generation  of  range-height  coverage  diagrams,  the  third  option  should  be 

selected,  by  entering  the  number  3  at  the  cursor  OPTION  prompt  Ohe  layout  of  the  first 

panel  is  shown  in  figure  4),  whereupon  the  panel  requesting  input  of  the  narameters  relevant 

to  the  generation  of  the  range  height  charts  will  be  displayed  (figure  5).  The  parameters 
are  self  explanatory,  with  the  exception  of  the  procedure  which  describes  the  atmospheric 
refractive  index  profile:  this  is  described  in  the  next  section.  Selection  of  GDDM  output 

will  display  the  range  height  chart  on  the  IBM  terminal,  and  selection  of  HP  output  will 

produce  Hewlett  Packard  graphics  language  (HPGL)  instructions  to  drive  a  Hewlett  Packard 

plotter  to  produce  a  hard  copy  of  the  chart,  on  the  output  file  ’JLW.RAYTRACE.PI  OT(name)’. 
These  instructions  may  be  subsequently  downloaded  to  a  HPGL  compatible  plotter  to  produce  a 
hard  copy  of  the  range  height  chart.  A  help  facility  is  available  via  the  PF1  key  (press 

once  for  short  message,  a  second  time  for  a  detailed  explanation)  if  the  user  is  uncertain 

of  the  appropriate  parameters  to  enter  in  any  of  the  panel  input  fields. 

3.2  The  procedure  specifying  the  refractive  index  profile 


The  RAYTRACE  package  contains  a  number  of  different  procedures  for  calculating  the 
refractive  index  of  the  atmosphere.  At  present  only  two,  namely  ATMOS1  and  ATMOS2  which 
contain  the  NBS  ngle  term  exponential  and  two  term  exponential  profiles  (ref.  5)  are 
appropriate  for  use  with  this  package.  These  profiles  are  respectively 


n  =  1  +  313.6  x  10  6e 


lh/7000) 


and 


,  .  n  ,,,-6  (h/9000)  ...  ...  6 

:i  =  1  +  265.0  x  10  e  +  48  x  10 

where  h  is  the  height  above  the  surface  of  the  Earth  in  metres, 
used  to  generate  the  range  height  cnart  examples  in  this  report. 


lh/2500) 

c 

The  procedure 


ATMOS  1 


was 


Any  alternative  refractive  index  profile  may  Ik  used  by  developing  another  procedure  to 
supply  the  refractive  index  as  a  function  of  height  The  call  to  the  procedure  is  of  the 
form 

CALL  ATMOS  (HEIGHT.  AI.PHAR,  REMI  X  SLOPE) 
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in  which  all  the  arguments  are  DECIMAL  FLOAT  (16)  variables,  (ie  double  precision).  The 
input  parameters  to  the  procedure  arc 

(1)  HEIGHT,  the  height  of  the  point  on  the  ray  above  the  surface  of  the  Earth 
(metres),  and 

(2)  ALPHAS,  the  angle  of  propagation  of  the  ray  with  respect  to  the  local 
horizontal  (radians). 

For  refractive  index  profiles  which  do  not  contain  any  discontinuities  (such  as  ATMOS  1  and 
ATMOS2)  the  parameter  ALPHAR  is  not  used,  but  is  included  in  the  CALL  ,o  provide 
compatibility  with  other  refractive  index  procedures.  The  procedure  then  returns 

(1)  REFIDX,  the  refractive  index  at  the  point  under  consideration,  and 

(2)  SLOPE,  the  rate  of  change  of  refractive  index  with  height  at  the  point 


Listing  1  gives  the  code  used  in  ATMOS  1,  and  may  be  used  as  an  example  to  generate  code  for 
any  other  refractive  index  profile. 

3.3  Examples  of  range  height-angle  charts 

A  lew  examples  of  the  charts  that  can  be  generated  by  this  software  are  included  to  show  the 
variation  in  shape  of  the  charts  with  different  values  of  the  input  parameters  that  might  be 

expected  when  using  this  software.  The  refractive  index  profile  described  in  ATMOS1  was 
used  for  all  these  charts.  Figure  6  shows  a  typical  chart  with  the  source  on  the  surface  of 
the  Earth,  and  the  maximum  range  and  height  such  that  the  ray  passing  through  the 
intersection  of  the  range  aid  height  boundaries  subtends  a  considerable  angle  at  the  source. 

The  second  chart,  figure  7,  shows  the  output  for  an  elevated  source  with  a  larger  range 
limit.  In  this  chart,  the  parameters  are  such  that  the  tangent  ray  to  the  Earth’s  surface 

lies  intersects  the  maximum  height  bourdary  before  the  maximum  range  is  reached.  As  a 
consequence  of  this,  the  chart  does  not  have  a  vertical  range  ''oundary.  Figure  8  shows  an 

example  with  a  larger  maximum  range  and  maximum  height,  in  which  the  tangent  ray  extends 

considerably  below  the  origin  of  the  chart.  Charts  such  as  this  one  provide  the  greatest 

difficulty  in  the  estimation  of  suitable  scale  factors  for  the  axes  so  that  the  resulting 
chart  completely  occupies  the  output  page.  The  final  example  (figure  9),  which  is  an 

expanded  portion  of  figure  8,  shows  the  increased  density  of  rays  (hat  the  software  includes 
near  die  zenith  when  the  scale  factors  are  appropriate. 

For  extreme  values  of  the  input  parameters,  in  particular  for  small  maximum  heights  in 
comparison  with  the  maximum  range,  the;  region  near  the  horizontal  ray  will  occupy  most  of 


i 
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[he  chart,  and  consequently  there  will  Ire  a  limited  number  of  points  fiom  which  to  construct 
the  range  and  height  curves.  In  such  ases  the  software  will  need  to  he  modified  to  provide 

a  greater  density  of  rays  between  say  2°  and  the  tangent  ray.  The  statements  requiring 

modification  are  in  the  DO  statement  controlling  the  angular  increments  (see  Listing  2) 
within  the  procedure  'JL’iV.RAYTRACE.PLI(MAIN3)'.  There  are  many  possible  ways  in  which  these 
statements  could  be  modified  depending  upon  the  specific  requirements.  Where  the  critical 
angles  are  known,  direct  replacement  of  appropriate  sections  of  the  IX)  specification  with 
actual  angular  limits  would  suffice,  and  would  be  the  easiest  to  implement. 

4,  CONCLUSIONS 

This  report  has  described  the  physical  principles  behind  a  computer  program  for  the 

generation  of  range-height -angle  charts  which  provide  an  indication  of  the  propagation  paths 
of  rays  launched  from  a  radar.  The  examples  giver  to  illustrate  the  use  of  the  program  have 
assumed  a  NBS  single  exponential  refractive  index  profile,  though  in  principle  this  program 
will  operate  with  any  refractive  index  profile  that  does  not  result  in  ducting  of  the 

electromagnetic  rays.  It  is  anticipated  that  minor  refinements  will  be  incorporated  in  the 
software  from  time  to  ,.me,  to  improve  the  qualify  of  output,  for  example,  but  no  further 
major  developments  are  pioposcd.  The  package  as  it  stands  is  a  useful  product  for  quick 
guidance  when  radar  coverage  questions  arise. 

ACKNOWLEDGEMENTS 


The  contribution  of  Mr  Tim  Miller  in  developing  the  user  interface  is  gratefully 
acknowledged. 


13 


SRL  0027-TM 


REFERENCES 

Ai.  Title 

Blata  1.  «  "Radio  Ray  (Radar)  Range- Height-Angle  Charts", 

NRL  Report  6650,  January  22,  1968 

SkolD'k,  "Introduction  To  Radar  Systems". 

Second  Edition  pages  446-7,  McGraw  Hill, New  York,  1980 

Baue".  "Range-Height-Angle  Charts  With  Lookdown  Capability". 

Microwave  Journal  vol  24,  no.  10,  October  1981 

Whitro  i.,.  and  "A  Radar  Propagation  Model". 

Miller,  3.G.  SRL  Technical  Memorandum  (to  be  published) 

Bean,  B.R.  and  "CRPL  Exponential  Reference  Atmosphere". 

Thayer,  G.D.  National  Bureau  of  Standards  Monograph  4,  October  1959 


SRL-0027-TM 


14  - 

ATMOS  1:  PROCEDURE  (HEIGHT, ALPHAR.REFIDX .SLOPE); 


/*  This  subroutine  is  used  to  determine  the  refractive  index,  and  */ 

/*  the  slope  of  the  refractive  index  at  the  specified  height  HEIGHT  */ 

/*  of  the  refractive  index  ptofile  specified  in  this  algorithm.  */ 

/*  The  parameter  ALPHAR  has  no  effect  except  at  a  discontinuity  of  */ 
/*  the  refractive  index.  ALPHAR  >  0  causes  routine  to  select  */ 

/*  refractive  index  approaching  HEIGHT  from  above,  <  0  from  below  */ 

/*  The  first  call  to  this  routine,  with  the  variable  INIT  set  equal  “/ 

/*  to  0  results  in  the  specification  the  constants  and  */ 

/*  discontinuities  in  refractive  index  profile.  Subsequent  calls  */ 

/*  will  have  INIT  set  equal  to  1  and  will  not  repeat  this  */ 

/*  initialisation  activity  Note  that  the  Earth’s  surface  is  treated  */ 

/*  as  a  discontinuity  in  this  software  so  there  wili  always  be  one  */ 

/*  discontinuity.  The  refractive  index  profile  specified  here  is  */ 

I*  the  standard  single  exponential  profile  with  no  additional  */ 

/*  discontinuities.  *7 


DECLARE  (HEIGHT, ALPHAR.REFIDX, SLOPE)  DEC  FL.OAT(!6); 
DECLARE  (NS.Z)  DEC  FLOAT(16)  STATIC, 

INIT  BIT(l)  INIT(’O’B)  STATIC; 

DECLARE  (DSCHGT(IOO)  DEC  FLOAT(16), 

NUMDSC  BIN  FIXED(3!))  EXTERNAL  STATIC; 

/*  test  INIT  to  determine  if  this  is  first  call  to  subroutine  */ 

IF  INIT  =  ’O’B  THEN 
DO, 

/*  specify  constants  for  calculating  refractive  index  */ 

NS  =  313. OE  6; 

Z  =  7000. 0E0; 

/*  specify  discontinuities  in  refractive  index  profile  */ 

NUMDSC  =  1; 

DSCHGT(l)  =  0.0E0; 

INIT  =  ’  1 '  B ; 

END; 

/*  calculate  refractive  index  */ 

REFIDX  =  1 .0E0  -t-NS*EXPf-HEIGHT'Z); 

/*  calculate  slope  of  refractive  index  V 

SLOPE  =  -NS*EXP(-HEIGHT/Z)/Z: 

END  ATMOS  1; 


Listing  1  Sample  Refractive  Index  Procedure  ( ATMOS  1) 


l 


DOANG  =  START  1  TO  END1  BY  INTV1, 

END1-1 .0  TO  MAXCUCANG,  10.00)  BY  -1.00, 
MAX(UCANG,  10.00)0.50  TO  MAX(UCANG,5.0)  BY  -0.50, 
MAX(UCANG,5.0)-0. 10  TO  0.00  BY  -0.10, 

-0.02  TO  TANANG  BY  -0.02, 

TANANG-0.1  TO  LCANG  BY  -0.10, 

LCANG-1 .0  TO  -10.0  BY  -1.00, 

START2  TO  END2  BY  INTV2; 


In  the  above  listing 

(1)  UCANG  corresponds  to  p  in  section  2.4 

(2)  TANANG  corresponds  to  fi  in  section  2.4 

(3)  LCANG  corresponds  to  y  in  section  2.4 


Listing  2 


DO  Loop  Parameter  Specification 
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Figure  3(a)  Chart  geometry  for  R  within  the  radar  horizon 

max 


Figure  3(b)  Chart  geometry  for  R  beyond  the  radar  horizon 
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-  - -  RAYTRACE  SELECTION  PANEL 

OPTION  -  =«  -=  > 

!  TRACE  RAYS  FROM  A  POINT  FOR  A  RANGE  OF  LAUNCH  ANGLES 

2  TRACE  RAYS  BETWEEN  TWO  POINTS 

3  PLOT  RANGE  HEIGHT  CHARTS 

4  PLOT  PROPAGATION  LOSS  VERSUS  HEIGHT 

USE  HELP  (PFI)  FOR  EXPLANATION  OF  INPUT.  PRESS  ENTER  TO  CONTINUE 
END  (PF3)  TO  TERMINATE  RAYTRACING  PACKAGE 


Figure  4  Raytrace  selection  menu 


RANGE  -  HEIGHT  CHART  PARAMETERS  - 


SOURCE  DATA 

ANTENNA  HEIGHT  =  =  =  >  (KFT) 

REGION  FOR  RAY  TRACING 

MAXIMUM  HEIGHT  ==  =  >  (KFD 

MAXIMUM  RANGE  =  =  ^>  (NM) 

RANGE  AND  HEIGHT  INTERVALS 

HEIGHT  INTERVAL  >  (KFT) 

RANGE  INTERVAL  =  =  =  >  (NM) 


SPECIFICATION  OF  ATMOSPHERIC  PROFILE  ROUTINE 
ATMOSPHERIC  PROFILE  =  =  -  > 

SPECIFICATION  OF  PROGRAM  OUTPUT  REQUIREMENTS 

GDDM  VDU  PLOT  =  =  =  >  (YES  OR  NO) 

HP  9872  PLOTTER  =  =  =  >  (YES  OR  NO) 

FILE  MEMBER  NAME  =  =  =  > 

USE  HELP  (PFI)  FOR  EXPLANATION  OF  INPUT,  PRESS  ENTER  TO  PROCESS  DATA 
END  (PF3)  TO  RETURN  TO  MENU 


Figure  5  Chart  parameter  specification  panel 
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